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Abstract. We investigate a model where neutrinos are strongly coupled to a new, 
light scalar field. In this model neutrinos annihilate as soon as they become non- 
relativistic in the early universe, and a non-zero neutrino mass has a marginal effect 
on the matter power spectrum. However, the angular power spectrum of the cosmic 
microwave background (CMB) is changed significantly because the strongly interacting 
fluid of neutrinos and scalars does not experience any anisotropic stress. Such models 
are strongly disfavoured by current observations. Interestingly, this leads to the 
conclusion that the relativistic energy density around the epoch of recombination must 
be in the form of very weakly interacting particles. This conclusion is independent of 
the specific interaction model. 
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1. Introduction 

The standard cosmological ACDM model with very light, weakly interacting neutrinos 
provides a very good fit to all current cosmological observations using only 6 parameters. 
Adding a significant neutrino mass of order 1 eV worsens the fit and in fact allows for 
strong cosmological bounds on the sum of all active neutrino rest masses. Within the 
ACDM model the upper bound at present is around J2 m u ^0.5 — 1 eV, depending on 
the specific priors and data used [1-11]. This limit is almost an order of magnitude 
stronger than the current upper bound on the effective electron neutrino mass from 
beta decays. Tritium decay measurements have been able to put an upper limit on the 
mass of 2.3 eV (95% conf.) [12], which translates into an upper bound on the sum of 
neutrino masses of roughly 7 eV. 

Very interestingly there is also a claim of direct detection of neutrinoless double 
beta decay in the Heidelberg- Moscow experiment [13, 14], corresponding to an effective 
neutrino mass in the 0.1 — 0.9 eV range. If this result is confirmed then it shows that 
neutrino masses are almost degenerate and well within reach of cosmological detection 
in the near future. 

However, the cosmological neutrino mass bounds do depend on a number of 
assumptions, and are thus not on quite the same footing as the direct laboratory 
measurements. Many different possibilities have been discussed for alleviating the 
cosmological mass bound. The conclusion is that the mass bound is surprisingly robust. 
For instance, including cosmological defects has little effect [15]. Another possibility, 
which has so far not been investigated, would be if the universe exits reheating at 
extremely low temperatures. In this case neutrinos would not fully thermalize before 
they decouple, and the end result would be that the neutrino number density is smaller 
than in the standard model [16-20]. 

Recently it was also proposed that neutrinos could be interacting strongly with 
a new, massless scalar field [21]. In this model, neutrinos annihilate as soon the 
temperature drops below the neutrino rest mass and transfer all their entropy to the 
scalars. The interesting feature of this model is that the neutrino rest mass never 
plays a significant role during structure formation so that cosmological neutrino mass 
bounds could disappear entirely. In the present paper we investigate this model in detail 
by solving the Boltzmann equation for a strongly interacting fluid of neutrinos and 
scalar. We find that, although the matter power spectra are indeed almost unaffected 
by the introduction of neutrino-scalar interactions, the spectrum of cosmic microwave 
background fluctuations is changed significantly. In the next section we discuss the 
formalism needed to study structure formation with a neutrino-scalar fluid. Section 3 
discusses a likelihood analysis based on present CMB and large scale structure data, 
and finally section 4 contains a discussion. 
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2. General scenario 



The model proposed in Ref. [21] introduces a new scalar field which interacts with 
neutrinos via a simple scalar or pseudoscalar coupling, 

C = hy Vi Vj (p + Qij Pij 5 Vj(f) . ( 1 ) 

Here we do not discuss the details of model building, but simply assume that such a 
coupling exists. In order to allow for almost complete annihilation of neutrinos at the 
present, the dimensionless coupling constant, g^, should be g^ > 10~ 5 for all neutrino 
flavours. 

On the other hand, as discussed in [21], couplings much larger than this are strongly 
disfavoured by neutrinoless double beta decay and by supernova observations. We 
therefore take g ~ 10~ 5 to be characteristic of the model. As was also calculated in 
Ref. [21], the annihilation rate of neutrinos in the non-relativistic regime is roughly 
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This should be compared with the Hubble expansion rate which is roughly given by 
H ~ 4.9 x 10 -12 h T e y 2 s _1 around recombination. 
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Figure 1. The neutrino annihilation rate at T = 0.3 eV for different neutrino masses 
for g = 10~ 5 . The horizontal dashed line is H at T = 0.3 eV. 

As can be seen, even for neutrino masses which are close to the upper bound from 
/3-decay measurements, the annihilation rate is much larger than the Hubble expansion 
rate until after recombination occurs at T ~ 0.3 eV. This means that in this model 
neutrinos and scalars can be treated as a perfect fluid with effectively zero mean free 
path for purposes of CMB and structure formation calculations. As will be discussed in 
the next section this has significant implications for CMB and structure formation. 
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3. The Boltzmann equation 

The evolution of any given particle species can be described via the Boltzmann equation. 
Our notation is identical to that of Ma and Bertschinger (MB) [22]. The simplest choice 
is to use synchronous gauge because the numerical routine for calculating matter and 
CMB power spectra, CMBFAST [23], is written in this gauge. As the time variable we 
use conformal time, defined as dr = dt/a(t), where a(t) is the scale factor. Also, as 
the momentum variable we shall use the comoving momentum qj = apj. We further 
parametrize qj as qj = qnj, where q is the magnitude of the comoving momentum and 
nj is a unit 3-vector specifying direction. 

The Boltzmann equation can generically be written as 

L[f] = ^ = C[f], (3) 

where L[f] is the Liouville operator. The collision operator on the right-hand side 
describes any possible collisional interactions. 

One can then write the distribution function as 

f{x\q,nj,T) =/ (g)[l + ^(^,g,n,,r)], (4) 

where fo(q) is the unperturbed distribution function. 

This applies to neutrinos as well as to the scalars, but with different / and *. 

In synchronous gauge the Boltzmann equation can be written as an evolution 
equation for * in /c-space [22] 
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where \i = n^kj and e = (q 2 + a 2 m 2 ) 1 ^ 2 . h and r\ are the metric perturbations, defined 
from the perturbed space-time metric in synchronous gauge [22] 

ds 2 = a 2 (r)[-dT 2 + (Sij + hi^dx^xP], (6) 
hij = J d 3 ke iU s (kikjhik, t) + (hkj - ^<y 6t7(£, r)) . (7) 

3.1. Collisionless Boltzmann equation 

At first we assume that j^C[f] = 0. This assumption will be relaxed later. The 
perturbation is then expanded as 

oo 

* = E(-0 , (2i + l)*ifl(|i). (8) 

1=0 

One can then write the collisionless Boltzmann equation as a moment hierarchy for the 
by performing the angular integration of L[f] 
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* 2 = jfe«-(2* 1 - 3* 3 ) -(-h+ -v) ^ (11) 
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It should be noted here that the first two hierarchy equations are directly related to the 
energy-momentum conservation equation. This can be seen in the following way. Let 
us define the density and pressure perturbations of the dark matter fluid as [22] 

5 = 5p/p (13) 

6 = ik 3 5T^/{p + P) (14) 

a ee -(^-^,)(T^-^T fe V3). (15) 
Then energy and momentum conservation implies that [22] 

^-(^HHgH* (16) 

e = £(1 _ 3w )e - -^e + - *V (17) 
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By integrating Eq. (9) over q 2 edq, Eq. (16) is derived and by integrating Eq. (10) 
equation over q 3 dq one retrieves Eq. (17). 

3.2. Collisional Boltzmann equation 

Starting from the general Boltzmann equation in synchronous gauge, Eq. (5), we now 
introduce interactions by lifting the restriction that j^C[f] = 0. Ideally, one should 
calculate the collision integrals in detail for some explicit interaction. However, for very 
strongly interacting neutrinos and scalars it is in fact quite simple to understand how 
the fluid behaves. 

The collision term in Eq. (9) is j dQj-C[f] and the one in Eq. (10) is j dQfij-C[f]. 
Integrating these two terms over momentum space one gets the collision terms in 
Eqs. (16-17) to be 



and 



C[f]dQq 2 dqe (18) 
J C[f]dQq 2 dqfiq = k l J C{f}dttq 2 dq qi (19) 

respectively. 

Looking at particle species i, for processes of type ii — > ii, any integral of the form 

J C[f]dttq 2 dqA, (20) 

where A e (I, e, is automatically zero because A is a collisional invariant. Thus, both 
the above integrals are zero, and the right hand side of the / = and 1 terms should be 
zero, reflecting that energy and momentum is conserved in each interaction. 
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All the higher order terms can be estimated from the relaxation time approximation 
[24], in which 

1 * 

T C[f] = --• (21) 

JO T 

Here, r is the mean time between collisions. Since the neutrino-scalar fluid is assumed 
to be very strongly interacting, the collision term dominates all other terms in the 
Boltzmann equation (since r <C H^ 1 ) unless \1> = 0. This argument applies to neutrinos 
as well as scalars, so that for I > 2, ty v j = = 0. In terms of the momentum- 
integrated equations, this means that a = 0, so that the fluid experiences no anisotropic 
stress. 

Since there are also creation and annihilation processes, ii — > jj, there will be 
non-zero terms for I = and 1, of the form [25] 

*«L = - *'•<>• < 22 > 

equivalent to the relaxation time approximation for ii — > ii, Eq. (21). Again, since 
t <C H^ 1 this means that these terms dominate unless = for 1 — 0,1. 

Altogether, this means that the combined neutrino-scalar fluid can be treated as a 
single fluid with a modified effective equation of state, 

w eS = ( 23 ) 

Pv + P<f> 

The perturbation equations to be solved are Eqs. (16-17) with w = w c s. This is exactly 
equivalent to the fluid of electrons, protons and neutral hydrogen around the epoch of 
recombination which can also be considered a single fluid because the interaction time 
scale is much shorter than the expansion timescale [34]. 

The effective equation of state is then straightforward to calculate. The temperature 
of the neutrino-scalar fluid, T V( ^, is found by solving Friedmann equation and the 
equation of energy conservation 

H 2 = 8nG(p u + p + p 7 + p m + p A ) ^ 

<n\ dL= _ H A P(j} + 3(p„ + p v ) 

alt al{ Pv + P4> )/alT v4 , ' 

with the energy density and pressure of each species given by 
3 f°° p 2 Ealp 



Pv 



± f°° Pbdp 
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P<t> = (28) 

P4, = (29) 

In Fig. 2 we show the effective equation of state of the combined neutrino-scalar 

fluid for a neutrino mass of 1 eV. At T > m F and T <^ m v the fluid is extremely 
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relativistic so that w — 1/3. During the annihilation epoch when the rest mass of the 
massive neutrinos plays a role, the pressure is lower. 
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Figure 2. The effective equation of state for the combined neutrino-scalar fluid for a 
neutrino mass of 1 eV. 

In Fig. 3 we show the energy density of neutrinos and scalars for the same model, 
normalized to that of a massless standard model neutrino. At high temperatures, 
p/Pu,o = 3 + 1 ~ 3.57, whereas after neutrino annihilation it is \ x | x (25/4) 4 ^ 3 ~ 6.58. 




Figure 3. The total energy density in the combined neutrino-scalar fluid for a neutrino 
mass of 1 eV. 
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4. Matter power spectra 

We have calculated matter power spectra for models with a fluid of strongly interacting 
neutrinos and scalars. 

Instead of plotting the power spectrum, P(k), we plot the transfer function T(k) 
which is defined according to P(k) = T 2 (k)P (k). Po(k) is the primordial power 
spectrum, which is independent of neutrino physics. We show these transfer function 
in Fig. 4 for a standard ACDM model with parameters Q = Q m + Q\ = 1, Q m = 0.3, 
Q b = 0.05, H = TOkms^Mpc- 1 , n s = 1, and r = 0. 

The stated mass is for a single mass eigenstate, and all mass eigenstates are assumed 
degenerate so that J2 m u — 3m u . Note that for m u = 0.01 eV this is unphysical since 
the mass difference 5m 2 3 is known from oscillation experiments to be of order 0.05 eV 
[35-38]. However, this is of no qualitative importance for the analysis. 
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Figure 4. Transfer functions for ACDM models with neutrino masses of 0.01 eV 
(solid), 0.1 cV (long-dashed), 1 eV (dashed), and 10 eV (dotted), respectively. 

From this figure it is clear that the matter power spectrum is not changed 
significantly, even for very large neutrino masses. This is in accordance with the finding 
in Ref. [21]. However, in Ref. [21] only the energy density in the neutrino-scalar fluid 
was accounted for. The strong neutrino-scalar interaction means that the fluid does not 
free-stream, but rather undergoes acoustic oscillations. 

In Fig. 5 we show the difference between the different power spectra on small 
scales. The difference on small scales between m v = and m u — > oo is roughly 
T(m u — > oo) /T(m u — 0) ~ 0.8, leading to the conclusion that the difference in the 
amplitudes of the matter power spectra is T 2 (m u — > oo) /T 2 (m u = 0) ~ 0.6. This is in 
almost perfect agreement with Ref. [21], where the ratio was found to be 0.65. 
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Figure 5. Transfer functions for ACDM models with neutrino masses of 0.01 eV 
(solid), 0.1 eV (long-dashed), 1 eV (dashed), and 10 eV (dotted), respectively. The 
normalization is arbitrary. 

5. CMB power spectra 

Contrary to the matter power spectra there is a significant effect on the CMB, depending 
on the neutrino mass. If the neutrino mass is large, neutrinos annihilate before matter- 
radiation equality, delaying it and thereby increasing the early integrated Sachs- Wolfe 
effect. 

However, there is also a significant effect from the fact that there is no stress 
term damping the neutrino-scalar oscillations. This effectively increases the fluctuation 
level on scales which are sub-horizon prior to recombination. In Fig. 6 this effect has 
been illustrated by simply setting all Z > 2 modes equal to zero for standard, massless 
neutrinos. Physically this corresponds to having neutrinos undergo acoustic oscillations 
instead of free-streaming. 

As can be seen from the figure, the effect is very similar to introducing a step-like 
feature in the power spectrum which would increase power on all scales smaller than 
^break, which should then be of order the horizon scale at decoupling. 

In Fig. 7 we show CMB spectra for same models as in Fig. 4. As can be seen 
there are several different effects interplaying. For the model with m v = 0.01 eV, 
recombination occurs long before annihilation, but for the models with m u = 0.1 and 1 
eV annihilation happens around the same epoch as recombination. Since the effective 
w is lower at this point there is a slight decrease in the early ISW effect compared to 
the m„ = 0.01 eV model. Finally, in the model with m v = 10 eV annihilation occurs 
long before recombination so that the early ISW effect is strongly enhanced because of 
the additional radiation present during recombination. 
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Figure 7. CMB TT power spectra for ACDM models with neutrino masses of 0.01 
eV (solid), 0.1 cV (long-dashed), 1 cV (dashed), and 10 eV (dotted), respectively. 
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The fact that the CMB spectrum is changed substantially makes it difficult to 
predict whether these models are in fact capable of fitting current data as well as the 
standard ACDM model. Since the CMB fluctuations are enhanced on sub-horizon scales, 
one prediction is that models with higher Q m and/or H will produce better fits since 
such models produce less power around the first acoustic peak. In the next section we 
present a likelihood analysis of present CMB and large scale structure data for models 
with a strongly interacting neutrino-scalar fluid. 

6. Likelihood analysis 

In order to test how well models with a strongly interacting neutrino-scalar fluid fit 
data, we have performed a likelihood analysis of CMB and LSS data. 

As our framework we choose the minimum standard model with 6 parameters: fl m , 
the matter density, the baryon density, H , the Hubble parameter, and r, the optical 
depth to reionization. The normalization of both CMB and LSS spectra are taken to be 
free and unrelated parameters. The only additional parameter is then the neutrino mass, 
m v . We assume that there are three active neutrino mass eigenstates with degenerate 
masses so that J2 m u = 3m u . 



Table 1. Priors on cosmological parameters used in the likelihood analysis. 
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Likelihoods are calculated from x 2 so that for 1 parameter estimates, 68% confidence 
regions are determined by A% 2 = x 2 — Xo — L an d 95% region by A% 2 = 4. Xo is X 2 
for the best fit model found. 

6.1. Large Scale Structure (LSS). 

At present there are two large galaxy surveys of comparable size, the Sloan Digital Sky 
Survey (SDSS) [28,29] and the 2dFGRS (2 degree Field Galaxy Redshift Survey) [27]. 
Once the SDSS is completed in 2005 it will be significantly larger and more accurate 
than the 2dFGRS. In the present analysis we use data from SDSS, but the results 
would be almost identical had we used 2dF data instead. In the data analysis we use 
only data points on scales larger than k = 0.15/i/Mpc in order to avoid problems with 
non-linearity. 
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6.2. Cosmic Microwave Background. 

The CMB temperature fluctuations are conveniently described in terms of the spherical 
harmonics power spectrum Cf T = (|a/ m | 2 ), where ^£-(9,(f)) = J2im a imXi m (9 , (j>) . Since 
Thomson scattering polarizes light, there are also power spectra coming from the 
polarization. The polarization can be divided into a curl-free (E) and a curl (B) 
component, yielding four independent power spectra: Cj T , C EE , C EB , and the T-E 
cross-correlation Cj E . 

The WMAP experiment has reported data only on Cf T and Cf E as described in 
Refs. [1,26,30-32]. We have performed our likelihood analysis using the prescription 
given by the WMAP collaboration [1, 30-32] which includes the correlation between 
different Cj's. Foreground contamination has already been subtracted from their 
published data. 

6.3. Results 

We consider two cases in the likelihood analysis: (a) Only WMAP temperature and 
polarization data, (b) WMAP data and SDSS data, combined with data on the Hubble 
parameter from the HST key project [33], h = H o /(100 kms" 1 Mpc" 1 ) = 0.72 ± 0.08. 
The priors we use are given in Table 1. 

In Fig. 8 we show x 2 as a function of m u , obtained from minimizing over all the 
other cosmological parameters. As can be seen from the figure the fit to data is quite 
poor in both cases when compared to the standard ACDM model with non-interacting 
neutrinos, independent of neutrino mass. Of course it might be possible to obtain a 
better fit by introducing additional parameters to the analysis, but this would be at the 
price of making the model more contrived. 

In Fig. 9 we show best fit values of Q m and h as functions of neutrino mass. As can 
be seen the best fit value of the Hubble parameter is highly non-standard for all values 
of the neutrino mass. In order to compensate for the increase in power due to the lack 
of free-streaming, the Hubble parameter has to be increased significantly in order to 
provide a decent fit. This in turn is in strong disagreement with measurements of the 
Hubble parameter from the HST key project [33] 

7. Discussion 

We have calculated in detail how structure formation proceeds in models with a 
strongly interacting plasma of neutrinos and massless scalars. These models differ from 
models with standard non-interacting neutrinos because neutrinos annihilate when the 
temperature drops below the rest mass. Furthermore, in these models the fluid of 
neutrinos and scalars undergoes acoustic oscillations instead of free-streaming. 

The effect is that the matter power spectrum is only minimally altered by a non- 
zero neutrino rest mass, but that the CMB power spectrum changes significantly, both 
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Figure 8. X 2 for WMAP TT and TE data only (full line), and WMAP+SDSS+HST 
data (dashed line). The horizontal lines show best fit x 2 for ACDM models with 
massless neutrinos. 
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9. Values of fl m and h for the best fit models, as a function of m„, for WMAP 
TE data only (full line), and WMAP+SDSS+HST data (dashed line). 
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because of the additional relativistic energy density and the strong coupling of the 
neutrino-scalar fluid. 

For standard, non-interacting neutrinos the neutrino free-streaming sets an upper 
limit on the neutrino rest mass when CMB and LSS data is combined. While neutrinos 
in the strongly interacting models do not suffer free-streaming and therefore produce 
matter power spectra which are in accordance with standard ACDM, it is impossible 
to simultaneously fit CMB and LSS data. Models with strongly interacting neutrinos 
and scalars are therefore strongly disfavoured by observations. Interestingly this also 
means that although present cosmological observations cannot distinguish between very 
light neutrinos and other types of relativistic energy density, it is possible to say that 
the relativistic energy density must be in the form of very weakly interacting particles 
which free-stream instead of oscillating acoustically (see also [39]). This is yet another 
confirmation of the standard picture of ACDM cosmology. 
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